function [i1,i2,T]=point(x0,we,mode,Te)
% global cur_d cur_q LdMap LqMap BAIDflux;
Rs=0.1;
eps=0.1;
k=1;
for i = 1:10
    [df,f] = fun(x0,we,mode,Te);
    x=x0-k*f/df;
    if(abs(x-x0) < eps)
        break;
    end
    x0 = x; % 更新迭代结果
end

i1=x0(1);
i2=x0(2);
% fai=interp2(cur_d,cur_q,BAIDflux,i1,abs(i2),'linear');
% Ld=interp2(cur_d,cur_q,LdMap,i1,abs(i2),'linear');
% Lq=interp2(cur_d,cur_q,LqMap,i1,abs(i2),'linear');
fai=0.038;
Ld=0.9e-4;
Lq=1.6e-4;
np=5;
T=1.5*np*(i2*(fai+(Ld-Lq)*i1));
